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We present an ab initio quantum theory of the metal-insulator transition in Ti2 03. The recently 
developed cluster LDA+DMFT scheme is applied to describe the many-body features of this com- 
pound. The conventional single site DMFT cannot reproduce a low temperature insulating phase for 
any reasonable values of the Coulomb interaction. We show that the non-local Coulomb interactions 
and the strong chemical bonding within Ti-Ti pair is the origin of the small gap insulating ground 
state of Ti2 03. 
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The complicated electronic structure and the nature 
of the metal-insulator transition in Ti2C>3 and V2O3 has 
been the object of intensive experimental and theoret- 
ical investigation over the past half century Re- 
cent progress in high-energy photo-emission spectroscopy 
and correlated electrons dynamical-mean field theory 
(DMFT) H has shed new light on the first-order metal- 
insulator transition in V2O3. It has been shown that an 
realistic description of the metallic and insulating phases 
of V2O3 can be obtained from the combination of a band 
structure scheme with the local electron-electron inter- 
action given from DMFT Q. The correlation effects in 
Ti 2 3 are less clear but angle resolved photo-emission 
experiment [f| shows a strong reduction of the Ti 3d- 
bandwidth compared to band structure calculations. The 
important question is related to the mechanism of the 
small, about 0.1 eV, semiconductor band-gap formation. 
The generally accepted view is that the metal-insulator 
transition is related to the decrease of the c/a ratio in 
rhombohcdral Ti203 and the formation of a Ti-Ti pair 
along z-axis Q . Below the broad (almost 250 K in width) 
metal-insulator transition at around 470 K the Ti-Ti pair 
distance is seen to decrease without any change of the 
rhombohcdral structure or the formation of long-range 
antiferromagnetic order Q • This is in contrast to the case 
of V2O3 where the V-V pair distance increases within a 
monoclinic distortion in the antiferromagnetic phase pj. 

Ti2C>3 has an CV-AI2O3 corundum structure (Fig. 
in the metallic and insulatingphases with two formula 
units per rhombohcdral cell |8j. Each Ti atom is sur- 
rounded by the octahedron of oxygens leading to the 
large t2 g -e g splitting. The trigonal distortion gives an ad- 
ditional splitting of tig bands into e g -ai g states and a\ g 
subbands of Ti-Ti pair form strong bonding-antibonding 
counterparts (Fig. In principle, the large decrease of 
the Ti-Ti distance could split further an occupied single- 
degenerate aig states from a double-degenerate e g states 
of tig subband and form the insulating d 1 configuration 
of this Ti compound. Nevertheless, state of the art LDA 
calculations have shown that for reasonable Ti-Ti pair 



distances Ti2C>3 will stay metallic 

In order to investigate the role of electron-electron 
interactions in the formation of this insulating low- 
temperature phase one needs an accurate estimation of 
the dig and e g bandwidths in this complex structure |lCj . 
For example a simple free [T^Og] 12- cluster mean-field 
investigation can easily produce a semiconducting gap 
due to drastic underestimation of the a\ g and e g band- 
widths O n the other hand a more accurate band 
structure calculation within the unrestricted Hartree- 
Fock approximation results in large gap antiferromag- 
netic state ^3 ■ Thus it is crucial to use both the correct 
Green-function embedding of the Ti-Ti pairs as well as 
a more accurate treatment of the electron-electron inter- 
action. 

The role of metal- metal pair formation and the " molec- 
ular" versus band pictures of the electronic structure 
have attracted much attention in these compounds |l3j . 
The combination of a strong on-site Coulomb interac- 
tion and the large anisotropy between the hopping pa- 
rameters in and perpendicular to the pair direction can 
favor a localized molecular-orbital picture of the insulat- 
ing phase. However, realistic tight-binding calculations 
for V2O3 have shown the importance of long-range hop- 
ping parameters ^J. It is also unclear how good an 
on-site approximation is for the electron-electron interac- 
tion. Since the pair forms a natural " molecular like" ele- 
ment in the corundum- type Ti203 structure it might be 
expected that non-local electron correlations are impor- 
tant in this system. Thus an approach which combines 
pair and beyond pair hopping with non-local electron in- 
teractions would be seem to be ideal for this problem. 

In this letter we apply for the first time a method, the 
cluster DMFT scheme 0, which contains all the 
physics of correlated pairs in crystals to determine the 
origin of the insulating phase and the metal-insulator 
transition in Ti2 03. A numerically exact multi-orbital 
Quantum Monte-Carlo (QMC) scheme is deployed for the 
solution of the cluster DMFT problem and an accurate 
first principles tight-binding parametrization used for the 
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one electronic structure. Our strategy here is to inves- 
tigate the gap formation using single site and clus- 
ter LDA+DMFT with only local correlations included. 
We then deploy the full non-local CDMFT and in this 
way are able to directly elucidate the impact non-local 
Coulomb interactions have on the physics. We show that 
the competition between strong bonding within the Ti-Ti 
pair and localization from correlation effects leads to the 
unique situation of the small semiconducting gap struc- 
ture in Ti203 oxide and that non-local Coulomb corre- 
lations are of crucial importance for the physics of this 
small gap insulators. 



LDA+DMFT scheme is calculated via the BZ-integration 
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FIG. 1: Left - rhombohedral unit cell of Ti2C>3 corundum 
structure. Titanium ions are indicated by the red color, oxy- 
gens by green, and the pair of Ti atoms in z direction by blue. 
Right - schematic representation of the tig splitting in Ti2C>3 
(top part) and the intersite Coulomb interaction matrix (bot- 
tom part). w=0 in all our calculations. 

We start with the orthogonal LDA Hamiltonian 
HH?i, (k) in the massively downfolded N-th order muffin- 
tin orbital representation |l8( {m corresponds to the 12 
tig orbitals of two Ti-Ti pairs in rhombohedral unit cell) 
and include different Coulomb interactions (see Fig. 
DMFT results for the local and non-local Coulomb inter- 
actions are presented in Figs. 12131 

The bare LDA density of states (DOS) is shown in 
Fig. El by the dashed lines for the low temperature struc- 
ture (LTS, ~300 K H) and high temperature structural 
(HTS, ~870 K ||) data on the upper and lower panels 
respectively. Both LTS and HTS electronic structures 
are metallic within the LDA scheme. The a\ g subband 
(green dashed line in the Fig. |3J has a strong bonding- 
antibonding splitting in contrast to the subbands (red 
dashed line). The bandwidth of the HTS is approxi- 
mately 2.8 eV and smaller than the bandwidth of the 
LTS (3.2 eV) due to the reduction of the t aigt a lg hopping 
from -0.85 to -0.63 eV. 

The cluster DMFT maps the many-body crystal sys- 
tem onto an effective self-consistent multi-orbital quan- 
tum impurity-cluster problem [T5L llrJ ] . The correspond- 
ing Green-function matrix for the Ti-Ti cluster in the 
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[(iuj n + fx)l - H LDA (k) - £K)]-\ (1) 



where [i is the chemical potential defined sclf-consistently 
through the total number of electrons. ui n — (2n + l)irT 
are the Matsubara frequencies for temperature T = /3 _1 
(n = 0, ±1,...) and a is the spin index. The Hamilto- 
nian and the self-energy matrix have the following super- 
matrix form corresponding to the symmetry of two Ti-Ti 
pairs in the unit cell 

( Hn + En H12 + S12 H13 H14 \ 

H21 + S21 H22 + Sn H23 H24 
H31 H32 H33 + Sn H34 + £12 

\ H41 H42 H43 + S 2 i H44 + En / 

(2) 

where Hy(k) and Sjj(w n ) are 3x3 matrices for the ti g 
states and En and E12 correspond to the intrasite and 
intersite contributions to the self-energy respectively. 

In the self-consistent cluster DMFT scheme the local 
Green-function should coincide with the correspond- 
ing solution of the effective two-site quantum impurity 
problem Q 

G a (r - t>) =-\\ D[c, c+]e- s ^c(T)c+(r') (3) 

here c(r) = [ci m(T (r)] is the super- vector of the Grassman 
variables, Z is the partition function, i runs over Ti-Ti 
pair and m runs over e g or a\ g orbitals. The effective 
cluster action S e f / is defined in terms of so-called "bath" 
Green function matrix Q C/~ 1 (w„) = G~ 1 (u; n ) + E CT (a;„) 
which describes the energy, orbitals, spin and tempera- 
ture dependent interactions of particular cluster with the 
rest of the crystal 

Seff = - / dr dr'Tr[c + {r)g- 1 (T,T , )c(r')}+ (4) 
Jo Jo 

2 E \Vmm' n i™ n im' + ^ mm' ~ " ^mm') n i™ n ■jm'] 



here rii mcT =c[ ma .Ci m(7 . We have parameterized the 
screened local Coulomb and exchange matrices (U£ m , 
and Jmm') f° r tne electrons in terms of average 
Coulomb and exchange integrals |l°J | and used a sim- 
ple approximation to the intersite , interactions as 
shown in the Fig. ^ 

The multi-band impurity QMC scheme jl^, 03] has 
been used for the numerically exact calculation of the 
cluster Green function (eq. . The number of auxiliary 
Ising fields in the discrete Hirsh-Fye transformation were 
15 and 58 for the local and non-local Coulomb interac- 
tion schemes respectively. For accurate QMC integration 
we used of the order of 10 6 sweeps, with 8000 k-points 
for the BZ-intcgration. Within 15-20 DMFT iterations 
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FIG. 2: Density of states for the single site (upper panel) 
and cluster (lower panel) DMFT calculations with different 
values of the Coulomb repulsion U and J— 0.5 eV. Inset: N(0) 
versus Coulomb parameter. Filled dark green circles - DMFT 
results with J— 0.5 eV, filled orange squares - DMFT with 
J—0 eV. Open blue circles and violet squares are CDMFT 
with J=0.5 and eV respectively. 

convergence in the self-energy was reached. The maxi- 
mum entropy method |22| has been used for analytical 
continuation of the diagonal part of the Green function 
matrix to the real energy axis. 

Firstly, in Fig|21 we show the total density of states 
for both conventional single site (DMFT) and the cluster 
(CDMFT) dynamical mean field theory where only lo- 
cal electron correlations have been included. The QMC 
simulation has been carried out for /3=20 eV -1 which cor- 
responds to a temperature of T~580 K which is on the 
border of the metal-insulator transition. In the upper 
panel of Fig. |2 are shown the results of DMFT calcu- 
lations with U =2,3,4 eV and exchange parameter J=0.5 
eV. For all values of Coulomb interactions there is a peak 
below the Fermi level at around -0.5 eV, predominantly of 
dig character with in all cases the same intensity. Above 
the Fermi level there are two peaks. The first is at 0.5 
eV and has e g character while the other peak is strongly 
dependent on the Coulomb parameter and can be asso- 
ciated with an upper Hubbard band. A lower Hubbard 
band can be seen at around -2 eV. We see that for all 
values of U the the shape of pseudogap is unchanged and 
the system remains metallic. On the lower part of Fig. |3 
the results of the CDMFT calculation are shown for the 
same values of the Coulomb and exchange parameters. 
The general structure of the DOS is seen to be similar 
to the single site calculation, however one may note in- 
teresting differences. The lower a\ g quasiparticlc band 
is decreased in intensity and shifted towards the Fermi 
level from -0.6 eV to -0.3 eV on increasing U from 2 to 4 



eV. This has the result that for U=4 eV the pseudogap 
is now located directly at Fermi level, whereas for other 
U- values and for all DMFT results it lies on the slope of 
quasiparticle peak. 

Using the temperature DOS at the Fermi level, de- 
fined as N(0) = ~ImG(uJo)/ir with ujo = irT we are able 
to estimate at what values of U the system will become 
insulating. This is indicated in the inset in the upper 
panel of the Fig. We see that for the single site cal- 
culations N(Q) depends weakly on U and the system will 
remain metallic up to very large values, about 8 eV, of the 
Coulomb parameter. On the other hand for the cluster 
calculation N(0) is seen to decrease strongly as a func- 
tion of U for both values of exchange parameter, and the 
critical value for an insulating solution is now lower at 
t/~5-6 eV. As expected for the d 1 configuration the fi- 
nite value of the exchange parameter effectively decreases 
the Coulomb interaction matrix. We see the single site 
results are in greater contradiction to the experiment as 
compared to LDA (see Fig. |2J : the local Coulomb inter- 
action leads to the reduction of the bonding-antibonding 
splitting of the ai g subband and this acts to suppress 
gap formation. On the other hand in the cluster case a 
small semiconducting gap is developed for large U due 
to dynamical antifcrromagnctic correlation in the Ti-Ti 
pair. 

Nevertheless, using cither the DMFT or CDMFT 
schemes with only local correlations there remains a dra- 
matic absence of gap formation in Ti2 03. We now deploy 
the full non-local correlation in CDMFT to the effect of 
non-local correlations on low and high temperature elec- 
tronic structure. We have used different values of the 
non-local Coulomb parameters and found that the most 
important correspond to non-diagonal interactions. For 
both structures we have chosen values of U=2 eV and 
J=0.5 eV which are close to those from constrained LDA 
estimations j2^| , while the off diagonal Coulomb param- 
eter W has been chosen at W=Q.h eV j^. On the upper 
panel of Fig. [3] is shown the total and partial DOS for 
/3=20 eV -1 . Shown also is the LDA result. One can see 
that for the reasonable parameters chosen we can repro- 
duce the correct value of the semiconducting gap ~0.1 
eV while keeping the bonding-antibonding splitting on 
the LDA level. In the lower panel the high temperature 
metallic solution corresponding to /3=10 cV _1 is shown. 
Here we emphasize that the proper inclusion of the struc- 
tural effect on the LDA level is important as evinced by 
the fact that for /3=20 eV -1 and high temperature hamil- 
tonian we again obtain a metallic solution. The e g states 
are similar for both LTS and HTS calculations with a 
small shift of occupied part in LTS case. However, the 
difference between the LTS and HTS phases is more pro- 
nounced for the dig states. The bonding-antibonding 
splitting in the LTS is about 2 eV whereas in the HTS 
case it is only 1.5 eV. The occupied a\ g states in the LTS 
phase are shifted down opening the insulating gap. The 
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FIG. 3: Partial and total CDMFT (solid line) compared to the 
LDA (dashed) DOS with W=0.5 eV and V aig =V eg =0. Total 
DOS are shown by black, the e g states by red, and ai g states 
by green. On the upper panel the low temperature struc- 
ture and /3=20 eV -1 are used. For the lower panel the high 
temperature structure and (3=10 eV -1 are used. The diago- 
nal and the biggest a\ g off-diagonal Green functions G(r) are 
shown in the upper inset. In the lower inset the _Re£ aisi a ls 
with W=0.5 eV are shown by blue, Re'E' aigjaig with W=0 eV 
by orange and ImE ai are shown by green 



important difference between the large U and small U 
plus non-local W is the absence of well defined Hubbard 
bands. This absence makes possible a critical test of the 
theory proposed here, and thus it would be very interest- 
ing for photo-emission experiments to check the existence 
or not of a lower Hubbard band at around -2 eV. 

We have shown that the cluster LDA+DMFT calcu- 
lation with a moderate Coulomb repulsion among the 
aig orbitals is essential to produce the high tempera- 
ture semimetallic state and the low temperature insu- 
lating state. To understand the role play of the inter- 
site Coulomb interaction we focus on the the quantity 
ta lg ,a lg + -ReE aig , ais (iw) which we can interpret as a fre- 
quency dependent "effective a\ g — a\ g hopping" which 
describes the hopping matrix element in the titanium 
pair. We find that this quantity is surprisingly frequency 
dependent (see lower inset of Fig. |3J) ■ 

We conclude that the main role of the intcrsite 



Coulomb interaction is dynamic (the Hartree contribu- 
tion to this quantity is small) and results in the effec- 
tive a\ g — ai g hopping that increases as the frequency 
decreases. This enhancement produces a strong level re- 
pulsion of the bonding antibonding ai g levels, lowering 
the ai g level relative to the e g level at the low frequency. 
This effect combined with a small narrowing of the ai g 
band opens the e g — ai g band gap which results in the 
insulating state. We checked that this enhancement of 
the effective hopping as frequency is decreased is absent 
if we turned off the intcrsite Coulomb repulsion. 

This effect is the cluster DMFT analog of a mecha- 
nism first discussed in the context of the single impurity 
model by Haldane [25J. He observed, that a Coulomb 
repulsion between the impurity level and additional con- 
duction electron states or screening channels, enhances 
the hybridization (single impurity analog of the hopping 
matrix element) as we renormalize to low frequency. 
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